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Abstract 

In the continuum the definitions of the covariant Dirac operator and of the gauge 
covariant derivative operator are tightly intertwined. We point out that the naive 
discretization of the gauge covariant derivative operator is related to the existence of 
local unitary operators which allow the definition of a natural lattice gauge covariant 
derivative. The associated lattice Dirac operator has all the properties of the clas- 
sical continuum Dirac operator, in particular antihermiticy and chiral invariance in 
the massless limit, but is of course non-local in accordance to the Nielsen-Ninomiya 
theorem. We show that this lattice Dirac operator coincides in the limit of an 
infinite lattice volume with the naive gauge covariant generalization of the SLAC 
derivative, but contains non-trivial boundary terms for finite-size lattices. Its nu- 
merical complexity compares pretty well on finite lattices with smeared lattice Dirac 
operators. 
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1 Introduction 



The standard mathematical description of the dynamics of the strong interactions, a 
description called Quantum ChromoDynamics (QCD), is obtained by writing the partition 
function of an Euclidean SU(3) quantum gauge field theory interacting with Nj > 2 
fermions in the fundamental representation of the gauge group. The partition function 
of QCD is postulated by analogy with the path integral formalism of Quantum Electro- 
Dynamics which has been shown to be successful with a very high accuracy. Using the 
rules of Grassmanian integration, the QCD partition function can be formally written as 
a functional integral over the non-abelian gauge degrees of freedom only, 




The measure of integration in eq. (1) can be interpreted as a formal probability measure 
over the space of gauge configurations because the Euclidean Dirac operator Df of each 
fermion flavor is antihermitian (with the right boundary conditions) and chirally invariant 
in the massless limit, 

D f = 7 A + m f > D n = d » + { 9K » ^ 
{7s, 7^} = » {7 M > 1v\ = 2^ , 7J = 7 M • 

Hence the eigenvalues of each operator Df come in complex conjugate pairs, up to a 
possible discrete set of zero modes of the operator 7^-D^, which guarantees reality and 
positiveness of their determinant for massive fermions, det.D/ > 0. 

The Euclidean gauge action (summation over repeated indices is implied throughout), 

S G = ^F^F^ , F^ = d^A v - d v A„ + ig [Ap, A v \ , (3) 

is invariant under the local gauge transformations G(x) G SU(3), 

A^x) — > G(x)A lx (x)G- 1 (x)+ l -(d,G(x))G-\x) , (4) 

whereas the operator Df transforms covariantly. 

Physical observables can then be related to expectation values with respect to the 
probability measure (1) of certain gauge-invariant matrix elements 0(A) of operators 
built out of the Dirac operators and their inverses, provided that the formal measure 
VA^ in (1) be given a precise gauge invariant meaning through a constructive procedure. 

Section 2 recalls briefly the main properties of Wilson's lattice regularization which is 
the only constructive proposal known to date [1]. The discretization of space-time in a 
finite box allows for the non-perturbative calculation of physical observables by means of 
numerical simulations. The approach has been very successful in describing many features 
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of hadronic physics, except for one thing. It proves difficult to reproduce the continuum 
physics with the physical light quark masses. 

The reason is well understood [3] and resides in the formulation of lattice fermions. 
The Lorentz invariant regularization of quantum fluctuations in continuum QCD generate 
a chiral anomaly which cannot be duplicated on the lattice with a discretization of the 
Dirac operator which is local, chirally symmetric and contains the correct number of 
fermionic degrees of freedom in the continuum limit. This result is usually referred to as 
the no-go theorem. The standard avoidance is to hold to local fermions and break chiral 
symmetry explicitly. 

In section 3 we reconsider the naive discretization of the covariant derivative and 
identify a set of local unitary operators which enables the definition of another lattice 
gauge covariant derivative with the same algebraic properties as in the continuum. This 
lattice derivative is non-local in accordance to the no-go theorem and coincides in the free 
limit with the SLAC derivative [4]. The lattice fermion formulation based on the SLAC 
derivative has been extensively discussed three decades ago. A general consensus has 
emerged according to which formulations of non-local lattice fermions coupled to gauge 
fields lead to various inconsistenties in weak coupling perturbation theory [5, 6, 7, 8, 9], 
and cannot reproduce the continuum limit properly, in particular the axial anomaly [10]. 
In fact we are not aware of a single numerical study of the functional integral of a SLAC- 
type fermion coupled to a compact lattice gauge field. 

Nonetheless we think it is fair to state that none of the objections to non-local fermions 
has the status of a no-go theorem. In all studies to date, the coupling of a SLAC fermion 
to a gauge field on the lattice has been written by mimicking the textbook derivation of 
a local gauge symmetry in the continuum. The result is correct only for infinite lattices. 
The boundary conditions on finite-size lattices are not taken into account by the standard 
technique. The unitary operators exhibited in section 3 are the right tools to include 
boundary conditions. As an example, in section 4 we diagonalize these unitary operators 
on periodic lattices and express in section 5 the matrix elements of the associated lattice 
Dirac operator in configuration space. We find non-trivial boundary contributions which 
vanish only in the limit of infinite physical volume. 

The constraints of the underlying locality of the exponentiated non-local lattice gauge 
covariant derivative probably cannot be neglected in analyzing the weak coupling pertur- 
bation theory, even in the infinite lattice volume limit. The intent of the present work is 
not to address this complicated issue, which deserves a separate work, but to stress its 
existence. 

In the concluding remarks we point out that the numerical complexity of the associated 
finite size lattice Dirac operator is similar to a five-dimensional local Dirac operator. 
In fact the algorithmic implementation is much simpler. Moreover this non-local Dirac 
operator can be naturally interpreted as smeared over the Wilson lines. This smearing has 
the virtue to be completely analytic. Numerical tests of the meaningfulness of the finite- 
size formulation of non-local fermions coupled to gauge fields can easily be performed on 
a single desktop computer in the quenched approximation up to four dimensions. 
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2 Lattice regularization 

As is well-known, Wilson's formulation consists of regularizing the Euclidean continuum 
gauge theory on a finite four- dimensional lattice C with hypercubic cells of spacing a, and 
of replacing the continuum gauge degrees of freedom, the gauge potential A^(x) which 
belongs to the SU (3) Lie algebra, by variables U x ^ associated to each link (x, x + ofi) of 
the lattice and which belong to the SU(3) group manifold. Then eq.(l) becomes 



where the integration measure is now a perfectly well-defined finite product of gauge- 
invariant Haar measures over the SU(3) group manifold, and Sg(U) and Mf(U) are 
discretized versions of the continuum gauge action and Dirac operators. This measure 
can be evaluated numerically by stochastic importance sampling. 

There is a large arbitrariness in the choice of lattice gauge action and lattice Dirac 
operators. The main constraint is that the lattice regularized model possess a second-order 
critical point which reproduces the asymptotic freedom of QCD in the continuum limit, 
with critical exponents predicted by perturbation theory. This requires in particular that 
the lattice operators reproduce the naive continuum definitions when the lattice spacing 
vanishes, a — > 0. Scaling theory then suggests that there exists a whole universality class of 
lattice actions which correspond to different regularizations of the same continuum theory 
and whose critical properties are related by renormalization group transformations. 

The most direct ab-initio approach is to consider lattice QCD actions with the same 
number of parameters as continuum QCD. The simplest such lattice gauge action, the 
Wilson action [1], has discretization errors of 0(a 2 ), 



Local gauge invariance is preserved on the lattice provided that the variables U x ^ trans- 
form as 



The simplest candidate for a lattice Dirac operator is expressed in terms of the naive 
discretization of the covariant derivative operator, namely, 






(6) 




(7) 



(8) 
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which has the correct covariant transformation law under (7) but does not have a spec- 
trum with definite transformation properties under the conjugation operation. Hence the 
determinant det T>i(U) is complex in general and does not define a probability measure. 

An obvious workaround would be to introduce antihermitian covariant difference op- 
erators, 

(D, t(l (U)i/>) x = — (u x ^ x+aP - Ul.tf^x-afi) , (9) 

which have the same conjugation properties as the continuum operators and produces 
a valid probability measure. But the operator j^D^ is plagued by the famous fermion 
doubling problem, due to the use of a central difference operator, and does not describe 
a single fermion flavor even in the continuum limit. 

Wilson proposed [2] to add to D s a piece proportional to the finite difference ap- 
proximation to the Laplacian operator A, which lifts the mass degeneracy of the fermion 
doublers by terms of order 1/a at the expense of breaking chiral invariance explicitly. 
However the Wilson operator still possess the same pseudo-hermiticity property as the 
continuum Dirac operator, 

V w (U)=^D s ^(U)-rA L (0<r<l), Vl = l5 V wlh , (fO) 

which guarantees the invariance of its spectrum under conjugacy and the interpretation 
of (5) as a probability measure. It was later realized [3] that it is not possible to devise a 
(ultra-)local lattice Dirac operator with the correct classical continuum limit and without 
fermion doublers while preserving exact chiral invariance on the lattice in the massless 
case. 

Various alternative lattice Dirac operators have been put forward during the subse- 
quent three decades, some of which with a presently viable ecosystem. The reader can 
find all references in a recent, and very nice, review [II] of the state-of-the-art of numerical 
simulations of lattice gauge theories. 

3 A chirally invariant lattice Dirac operator 

If one examines the definition (8) of the naive lattice covariant derivative operator, one 
realizes immediately that the four operators, 

Sfj, = 1 + a Di tfJ/ , (SijJixjy — (U x ,x+ati)ij ^y,x+ap , (11) 

are unitary ladder operators which translate by one lattice unit in direction fi each slice 
of the lattice field they act upon, while rotating locally their color degrees of freedom, 

fl^ = = 1 , V/x. (12) 

The set of operators transforms covariantly under the local gauge transformations (7), 

S„ — > GS„g-\ (0) ix j y = G(x) i:j Sxy , V/x, (13) 
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and encodes all the space-time and color degrees of freedom of a gauge field configuration 
on a four-dimensional lattice. For instance, the Wilson action (6) can be written, up to a 
constant term, as 

S W (U) = -^Tr {S^SlSt) . (14) 

Expressing the unitary operators S^U) as exponentials of antihermitian operators D r ^(U), 

S,(U) = e aD ^ u \ D r ,, + D\„ = , V/i, (15) 

singles out the operators D r ^(U) as the natural definition of the lattice covariant deriva- 
tive. Indeed the symmetric covariant difference operators (9) are just the leading approx- 
imation in the series expansion of these exponentials with respect to the lattice spacing 

a, 

aD s jU) = l - (e aD - - e- aD -) . (16) 

Then we can define the lattice Dirac operator 

V r (U)=^D r ^(U)+m, (17) 

which is antihermitian, chirally symmetric in the massless limit, and transforms covari- 
antly under the local gauge transformations (7). The eigenvalues of the operator T> r come 
in complex conjugate pairs, m ± i\, up to a possible set of zero modes for the imaginary 
part which ensures, like in the continuum, reality and positiveness of the determinant for 
massive fermions, detV r > 0. 

The operator V r (U) is non-local since each lattice covariant derivative operator D r ^ 
is a series expansion in the local covariant derivative operator 

aD r . tll = log(l + aD hjl ) 

- 2^ - v i# ■ 

71=1 

From the convergence properties of the series expansion (18), the prospect of a practical 
numerical implementation of the lattice Dirac operator T> r might seem very slim. 

On the other hand, in the free case, the eigenvectors of the operators are just plane 
waves and their 3iV 3 degenerate spectrum is simply 

Spec 5,(1) = je-" , p, = ^ , < n < ^} . (19) 

In this limit the lattice Dirac operator T> r has a discrete Fourier representation which has 
the same form as in the continuum, 

Vr(l) = ilvP» + m, (20) 
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which shows that the operator V r is not afflicted with the fermion doubling problem. In 
compliance to the Nielsen-Ninomya theorem, the price to pay is the non-locality of the 
operator T> r (U). The operator D r ^{\) coincides with the SLAC derivative introduced 
long ago [4] which reads, in the limit of an infinite lattice volume, 

Do^{x-y)= [*' a **i Plt e**-yK (21) 

J-TT/a {^) 

The usual recipe to couple a SLAC-type fermion to a SU (3) gauge field consists in restor- 
ing gauge invariance by inserting, if x — y has a non- vanishing component only in direction 
/2, the Wilson SU(3) ordered straight line integral W^(U,x,y) between x and y, 

(Ds(U)) MxMy = m 5 aP 5 l3 5 xy + £ ( 7 „)^ D^(x - y) (W,(U, x, y)) tj , (22) 

(y M -z M -l)/a 

W^U, x, y) = <5 X ± y ± ] ] U x+kan , x+{k+1)an , if y^ > x^ . (23) 

fc=0 

where x = (x^, x^), y = (y^, y^), and x^, label the sites in the three-dimensional slices 
orthogonal to the /ith direction. If y^ < x^, we have of course W^(U, x, y) = W^U, y, x)K 

The matrix elements (22) are certainly correct in the limit of an infinite lattice. How- 
ever the boundary conditions on finite-size lattices are not taken into account by the 
conventional prescription. The operators have such a simple structure that it makes 
possible their explicit diagonalization for an arbitrary background lattice gauge field con- 
figuration and arbitrary boundary conditions. We shall find non-trivial boundary terms 
in the matrix elements of the operator T> T on finite-size lattices. 

To the best of our knowledge, the underlying local unitary structure of the lattice gauge 
covariant generalization of the SLAC derivative does not seem to have been appreciated 
since its introduction. 



4 Explicit diagonalization of the unitary operators 

For definiteness we shall assume the lattice to be hypercubic, N\ = N 2 = N 3 = N 4 = N, 
and we shall impose periodic boundary conditions on both the lattice gauge field config- 
uration and the lattice matter fields which can be scalars, fermions, • • • . The operators 
do not commute in general and have different eigenspectra {A^, ifix}. Iterating the 
eigenvalue equation for S^, 

S^x, = A^a, , (24) 

yields 



xk=0 



j,x+nan ' £25) 
ij 
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Imposing the periodic boundary condition, 



ipx+Nap = ip x , Vx , (26) 
implies that the eigenvectors ^ satisfy the equations 

(W, (U, x,x + Nap)\ 3 {^). x = X" (^) ix , V* , (27) 

where (U, x,x + Nap,), defined in (23), is the Wilson line from x in direction p which 
wraps the lattice. We shall use the shorthand W^ X {U) for such Wilson lines and the 
argument U will be implicit most of the time. 

Hence each non-zero space-time component of the eigenvector ^ is a color triplet 
which is an eigenvector of some Wilson line W^ X (U). The Wilson lines are covariant 
objects under the local gauge transformations (7) and their eigenvalues are gauge- invariant 
and do not depend on the choice of base points x which differ only by the x^ component 
along their direction. Indeed a change of base point along a Wilson line is nothing but a 
similarity transformation. 

Therefore the eigenspectra of the Wilson lines W^ X (U) can be labelled by the points 
x of the lattice slice x M = 0, x = (x^, 0), 

{e^ , rfe} , -7T < 5l s < vr , c = 1, 2, 3 . (28) 

where rf s are the color triplet eigenvectors of W^(U). Their calculation requires only 
4iV 4 £77(3) matrix multiplications and 4iV 3 SU (3) matrix diagonalizations. 

Barring accidental degeneracies, the eigenvalues A M and eigenvectors fall into fami- 
lies labelled by the 4 x 3 x N 3 eigenvalues of the Wilson lines and defined by the equations, 

A^ = e iS °^ . (29) 

The general solution for A M reads 

K^-*******"", P„-^-£, 0<»<AT. (30) 

The N non- vanishing components of the corresponding eigenvector i/)° Sp are, with n = 
().••• .A - 1. 

with Vrf+vpji = w t (x, x + y^p) rf^ . 

These components can be computed sequentially and the calculation of each of the 127V 4 
eigenvectors requires only N £77(3) matrix-vector multiplications. So the total computa- 
tional complexity of the complete diagonalization of every operator £^ is of order 0(N 5 ). 
A complete diagonalization has to be performed only once for each lattice gauge field 
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configuration and its storage requirement is in practice proportional to the lattice volume 
since it is more efficient to recompute the N components of each eigenvector when needed. 

The eigenvectors of are also eigenvectors of D r ^ and the action of D r>II on an 
eigenvector 1])°^%^ has a remarkably simple continuum-like expression, 

aD r,u%,x,p» = i { a Pn + %,s, P ^ , = ■ ( 32 ) 

There is an inherent ambiguity in the logarithmic definition of D rfl . We have defined 
rather arbitrarily the phases 8^ - as the principal argument of the eigenvalues of Wilson 
lines. Other prescriptions are possible for unitary gauge groups. 



5 Matrix elements of the operator V r (U) 

The kernel operation which enters most algorithms involving fermions, such as the calcu- 
lation of the fermion propagator, is the action of the lattice Dirac operator on an arbitrary 
lattice fermion field. The action of T> r on a fermion \& can be written spin-component-wise 

as 

(2> r (E0*) a = (7m) q/3 DrAU)*? + m^ a . (33) 

To perform this calculation we just need to expand each spin component of the fermion 
field over the complete eigensystem of every operator D r4l , 

^ = E C l^,v.%^ > V ^ • ( 34 ) 

C,X,Pn 

We get 4 x 4 x 3 x N 4 equations, with x = (x, x^) and a = 1 throughout this section, 

(?P)jx = E^./Wm ( A ^,Pm) WU)j ■ (35) 

We can always choose all color triplet eigensystems to be orthonormal. Then 

we observe that all eigensystems {r)° x } along the same Wilson line are simultaneously 
orthonormal, 

Eto),«4-=^' Vx = (x,x /i ), (36) 
3 

since they are related by unitary transformations which preserve the scalar product. Thus 
we can transform each equation (35) into a simple one-dimensional Fourier series, 

W)^e-^^te).(^, 

^ 3 (37) 
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Therefore the coefficients Cp fijll ^p are the one-dimensional inverse discrete Fourier trans- 
forms, 

(38) 

It follows that the total computational complexity of the action of the operator V r (U) on 
a fermion field is of order 0(N 5 ) which is only a factor N more expensive than the action 
of a a local operator like the Wilson operator T> W (U). 
Plugging (38) into (34) yields, with x = x + x^fi, 

( D rM iV =jf E ^+%^- l(P ^ sW ( 39 ) 



c,x,x 



Inserting (31) gives, with y = y + y^fi, 

(40) 

The summation over p M brings in the finite-size SLAC derivative, 

^w = ^E^ e ^. ( 41 ) 



N 



and the summation over x produces, 



= s^dnAv,-^) E (E e ^ ( ^ M) to), to) J (*/») 
+*E [ to) J <*f>)ky 



kx 

(42) 



fc \ c 

The last step is to insert (31) and use the identity, 

(Wj (y, y + y^)) jt (E f (& 9 ) «v)i fe) m ) W (£ V + ifcflU = (/ (W*v))* ■ 

(43) 

Collecting everything together, the matrix elements of the operator V r (U) read finally, 

V m / (44) 

x E (^)^^(y, - *„) (w,(U,y,x) (W,, x (U))^- x ^ Na ) , 
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where we reintroduced the lattice spacing a to make apparent the physical dimensions. 
The covariance of the operator V r (U) under the local gauge transformations (7) is clearly- 
satisfied, 

v r (u) — >■ gv r (u)g- 1 . (45) 

We find two additional boundary contributions with respect to the infinite volume ex- 
pression (22). The first one is a diagonal term in configuration space which vanishes 
proportionally to the inverse physical lattice size. The second one is an insertion in the 
open Wilson line Wp(U, y, x) of the closed Wilson line W^ X (U) raised to a power the vari- 
ation of which is also proportional to the distance y — x in physical units. The insertion 
point can be covariantly transported anywhere along the closed Wilson line. 

6 Outlook 

As recalled in the introductory section, the lattice fermion formulation based on the SLAC 
derivative has been rather controversial. There have been one-loop calculations with the 
SLAC operator (22) in weak coupling perturbation theory of lattice Quantum ElectroDy- 
namics (QED) in four dimensions, which have shown the occurence of singularities in the 
fermion triangle graph [5] and in the vacuum polarization [6], that lead to non-local, non 
Lorentz covariant expressions. These singularities are generated by the discontinuities of 
the SLAC derivative at the edges of the Brillouin zone, = ±ir/a. It has been claimed 
[7] that these divergences could not be renormalized while keeping in the continuum limit 
a — > (at L = Na — > oo), both chiral invariance without extra states and Lorentz 
invariance. However it has also been suggested [8] that QED could be recovered in the 
continuum limit by a proper, non-perturbative, treatment of the infrared singularities and 
by imposing a finite number of non-local renormalization conditions. But the application 
of such empirical prescriptions to the lattice Schwinger model, namely two-dimensional 
Quantum Electrodynamics with massless fermions, which is a completely solvable model 
in the continuum [12], has still generated a spectrum doubling, a vanishing anomaly, a 
vanishing vacuum expectation value for (tpip), and a non-covariant axial- vector [9]. 

Despite all these negative results, a rigorous treatment of the perturbative expansion 
of the lattice Dirac operator (22) in the infinite volume limit is still lacking. As already 
emphasized, the underlying unitarity structure of the non-local lattice covariant derivative 
has not been taken into account in existing studies and the issue of spectrum doubling 
in such a formulation, which depends crucially upon the handling of singularities in the 
infinite volume limit, has to be settled accordingly. 

Whatever the outcome, we advocate a pragmatic approach, a-la Wilson. We have 
exhibited a covariant and chir ally- invariant lattice Dirac operator on finite-size lattices 
which has certainly no spectrum doubling in the free limit. Moreover, the expressions (44) 
of the matrix elements of the operator V r are quite convenient for an actual computer 
implementation. If spectrum doubling does occur when the gauge interaction is turned on, 
one could always add an explicit chiral symmetry breaking term. Our derivation reveals 
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that the non-local operator V r is a smeared operator with a controllable analytic averaging 
over the links of Wilson lines (whereas the original smearing proposal [13], as well as its 
many variants, are empirical thickenings of the links). Smeared operators are smoother 
and it is widely known [11] that their inversion has better convergence properties than 
local operators. So the accelerated convergence near the chiral limit, whose qualitative 
nature can be studied in the quenched approximation, may even turn out to compensate 
for the additional computational complexity of V r with respect to local Dirac operators 
such as V w . 
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